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Abstract 

We compute the hard photon production rate of a chemically lion-equilibrated 
quark-gluon plasma. We assume that the plasma is already thermally equili¬ 
brated, i. e. describable by a temperature, but with a phase-space distribution 
that deviates from the Fermi/Bose distribution by a time dependent factor 
(fugacity). The photon spectrum is obtained by integrating the photon rate 
over the space-time evolution of the quark-gluon plasma. Some consequences 
for ultrarelativistic heavy ion collisions are discussed. 

I. INTRODUCTION 

Hard photons are a promising probe for the fireball created in ultrarelativistic heavy ion 
collisions [1], As they leave the fireball without further interaction [2] they are probing the 
various stages of the collision directly. Concerning the quark-gluon plasma (QGP) phase, 
photon emission has been considered so far mostly at thermal and chemical equilibrium [3-5]. 
Here we will investigate the photon production from a chemically non-equilibrated QGP, 
starting from a simple model for the chemical equilibration of a parton gas in ultrarelativistic 
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heavy collisions [6]. This model predicts that the QGP, possibly formed at RHIC and LHC, 
is always far from chemical equilibrium, showing a strong undersaturation of the phase space 
in particular for quarks, which are the source of photon emission. In the next section we 
present a short review of the parton chemistry model, before we calculate the hard photon 
production rate from a chemically non-equilibrated parton gas in section 3. In section 4 we 
show photon spectra, obtained by integrating the photon rate over the space-time evolution 
of the fireball, and discuss some consequences for RHIC and LHC. Similar investigations 
have been performed in ref. [7-9], to which we compare our results. 


II. PARTON CHEMISTRY 


The chemical equilibration of quarks and gluons has been described by means of rate 
equations for the quark and gluon phase space density after thermal equilibrium set in [6]. 
Here we speak of thermal equilibrium as soon as the momentum distribution of the par- 
tons becomes exponential and isotropic. According to the event generator HIJING [10] the 
primary hard parton collisions result already in an exponential ^-distribution [11], Sub¬ 
sequent longitudinal expansion leads to an isotropic momentum distribution in the central 
slice (Az = 0.5 fm) at a time r 8S0 = 0.3 fm/c at RHIC and r 8S0 = 0.2 fm/c at LHC after 
the primary collisions [6] corresponding to r 0 = 0.5 — 0.7 fm/c after the maximum overlap 
of the nuclei [11,12], 

For times r > r 8S0 we assume that the distribution functions can be approximated by 


n F,B(E ) \,gi T ) e E/T(T) z j_ q ’ 


( 1 ) 


i.e. by equilibrium Fermi/Bose distribution functions with a time dependent temperature 
multiplied by a time dependent factor \ q ^ g which describes the deviation from chemical 
equilibrium. This factor, called fugacity, takes account of the undersaturation of the parton 
phase space density, i.e. 0 < \ q ^ g < 1. 

To lowest order perturbative QCD the phase space will be populated by the reactions 
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gg ^ ggg, 


gg qq- 


( 2 ) 


The time evolution of the quark and gluon densities, which are proportional to the fugacities, 
can de determined by rate equations, where the equilibration rates entering these equations 
follow from the cross section of the above reactions. The cross sections are calculated from 
the lowest order matrix elements, where the thermal gluon and quark mass, also depending 
on the fugacities, serve as infrared cut-offs: 
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( 3 ) 

( 4 ) 


where Nj denotes the number of active flavors in the parton gas, for which Nj = 2.5 was 
chosen. (In chemical equilibrium the thermal masses follow from the zero momentum limit 
of the gluon and quark self energies in the high temperature approximation [13] leading to 



m «=i£ r dkk{ ’ '3(‘) +n ’(‘)]' 


In the chemically non-equilibrated parton gas the equilibrium distribution functions are 
replaced by (1) leading to (3) and (4).) In addition, gluon radiation in the first reaction is 
assumed to be suppressed if the mean free path of the gluons in the parton gas is smaller than 
the formation length of the emitted gluon. In this way the so called Landau-Pomeranchuk 
effect is taken into account phenomenologically. The results for the equilibration rates can 
be found in ref. [6]. 

Together with the equation describing energy conservation in the case of a purely lon¬ 
gitudinal expansion [6] the rate equations determine the evolution of the quark and gluon 
fugacities A g , 3 (r) and of the temperature T(r). These equations are solved numerically 
together with the following initial conditions at r 8S0 resulting from HIJING: A° = 0.09, 
A° = 0.02, and T 0 = 570 MeV for RHIC and A° = 0.14, A° = 0.03, and T 0 = 830 MeV for 
LHC, respectively. The results are shown in fig. 1. (Recent investigations [12] with slightly 
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changed initial conditions and a somewhat larger rate for the gluon production gave similar 
results). The temperature drops faster than in the Bjorken scenario [14] (T 3 r=const), since 
energy is consumed by parton production. The main result of these investigations is a clear 
deviation from chemical equilibrium (undersaturation) at RHIC as well as LHC, especially 
for quarks. 

The following problems and criticisms associated with this model and its results should 
be mentioned: 

1. The initial fugacities from HIJING are very small in contrast to the one following 
from the Parton Cascade Model [15], which gives A° ~ 1 and A° ~ 0.7. The reason for this 
essential difference is not clear yet. 

2. The introduction of time dependent distribution functions such as (1) for describing 
a non-equilibrium situation contradicts a perturbative expansion in the real time formalism 
by leading to singularities in loop diagrams [16]. However, here we only want to employ the 
distributions (1) as a phenomenological ansatz, and we do not consider loop diagrams. 

3. Uncertainties in the equilibration rates have a significant influence on the evolution 
of the fugacities. In particular, a larger gluon production rate could be obtained if the 
Landau-Pomeranchuk effect in the gluon bremsstrahlung is treated by taking into account 
the rescattering of the radiated gluon. Then gluons with larger formation times can also be 
emitted [17]. Furthermore, the process gg ^ n g 7 n > 3, can be included leading to gluon 
production rates about twice as large [18]. 

Despite all these uncertainties, we think it worthwhile to study the consequences of 
this parton chemistry scenario, as it has been done in the case of charm production [12], 
The other way round, by comparing the particle production (photons, dileptons, charm), 
predicted by this scenario, with experimental data, one might be able to extract information 
on the equilibration in ultrarelativistic heavy ion collisions. 
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III. PHOTON PRODUCTION RATE 


In thermal and chemical equilibrium, the production rate of hard photons with energy 
E T can be computed using the Braaten-Yuan prescription [19]. This results in a decom¬ 
position of the rate into a soft part, which is treated using a resummed quark propagator 
according to the Braaten-Pisarski method [20], and a hard part containing only bare propa¬ 
gators and vertices. In the soft part the resummed quark propagator takes care of medium 
effects in the QGP; e.g., it contains the thermal quark mass m 2 = g 2 T 2 / 6, which serves 
as an infrared cut-off in the case of a vanishing bare quark mass. The hard part follows 
from the momentum integration over the matrix elements that lead to photon emission in 
lowest order (quark pair annihilation, Compton scattering with an initial gluon) multiplied 
by the distribution functions of the incoming and outgoing partons [3]. A separation pa¬ 
rameter k c is introduced, which allows to distinguish between soft and hard momenta of the 
intermediate quark. Assuming the weak coupling limit and gT <C k c <C T , the final result 
is independent of the separation scale k c . This procedure has been demonstrated in refs. 
[3,4] using Boltzmann distribution functions for the incoming particles, and in ref. [5] using 
full Fermi/Bose distribution functions with a nonzero quark chemical potential g. Unfor¬ 
tunately, the Braaten-Pisarski method is based on the principle of detailed balance, which 
holds only in full equilibrium, and thus is not applicable for the chemical non-equilibrium 
stage. 

Since there exists so far no consistent method for treating medium effects at non¬ 
equilibrium, we propose the following procedure. We only consider the hard part of the 
photon rate and replace the cut-off k 2 by 2m 2 as suggested by the equilibrium result [3]. 
For the thermal quark mass we adopt the formula (4) containing the fugacities, thus taking 
into account non-equilibrium effects. This approximation is in line with the estimates of the 
equilibration rates calculated in ref. [6]. In addition, we assume g = 0 since the photon rate 
is not sensitive to a non-zero quark chemical potential [5]. Furthermore, the fugacities show 
up in the distribution functions of the incoming and outgoing partons, for which we adopt 
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the non-equilibrium distributions (1), under the momentum integral defining the photon 


We now compute the hard photon production rate, starting from the following equation 
for each contribution (annihilation, Compton) of the photon rate [5] (summed over the 
photon polarisations): 

2E i^r= 8?2-W 7 * / <u £ i* i 2 / iE ■ ^ 9(F(£i f; )>ni,,2 | 1±n3> . (s) 

il “ ,X mw)E 'M -U |R> JP(EuE>) 

where )T) \Ai \ 2 denotes the square of the matrix element for the annihilation or Compton 
process summed over the initial and final parton states and rii the parton distributions, 
where the plus sign in front of n 3 corresponds to the annihilation process and the minus sign 
to Compton scattering. The polynomial P is given by P(Ei 7 E 2 ) = — {tE\ + (s + t)E 2 ) 2 + 
2 Es((s + t)E 2 — tEi) — s 2 E 2 + s 2 t + st 2 with the Mandelstam variables s, t, and 0 is the step 
function. E and p are the outgoing photon energy and momentum, related by E = p = |p|, 
since the photon is on-shell. Finally, dn is the number of photons that will come out of the 
plasma cell d 3 x during dt with a momentum in d 3 p around p. Equation (5) is already written 
in a Lorentz invariant fashion (for boosts along the beam direction) and can be related to 
the rapidity via E/d 3 p = dp x dp y dy = d 2 pjdy. 

Let us now introduce the fugacity factors in (5) by the replacement 

n 1 n 2 (l±n 3 ) X 1 n 1 X 2 n 2 (l ± A 3 n 3 ) . 

We decompose this product as follows: 


AiniA 2 n 2 (l ± A 3 n 3 ) = X 1 X 2 X 3 n 1 n 2 (l ± n 3 ) + X 1 X 2 (1 - X 3 )n 1 n 2 . ( 6 ) 

Here, the first term leads to the equilibrium photon rate as computed in ref. [3] multiplied 
by A^ X g \ 

( dn \ 'oaa s X 2 X q 9 W/T T ( 4ET\ . . 

( 2E ——— I = -1_jl^ T 2 e - E /T ln —— -L42 . 7 

y d 3 pd 4 x J x 97 t 2 |_ y k 3 J \ 

Here k 2 is the infrared cut-off of the hard contribution for which we will use 2 m 2 as dis- 

c q 

cussed above. In order to find (7), one has to use Boltzmann distribution functions instead 
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of the full quantum-mechanical distribution functions on the incoming legs. This underesti¬ 
mates the contribution coming from the Compton scattering process and overestimates the 
contribution of the quark-antiquark-annihilation process, but the two errors cancel up to a 
maximum error of about 10% in the end, as is shown by numerical analysis in ref. [5]. 

The second term of (6) results in the rate 

/ oj? ^ n ^ ^1^2(1 ~ ^3) 

\ d 3 pd 4 x J 2 8(27 t) 7 E 

comp 

,j d s J it J i. (8) 

2kl -s+kl E+S/4E |R V r \ £j +^ rj ‘ 2 > 

We treat this term again in the Boltzmann approximation, using nin 2 = e~ E+ ^ T . The 


distribution functions are constant in the innermost integral, and can be drawn out in front. 
Now all the integrals are elementary, and (8) reduces in a few steps to 



dn \ 
d 3 pd 4 xJ 


[A»A s (l-A,)[l+2eW''“ T E 1 (J% 


+A 9 A,(1-A i7 ) -2 + 2e 11 7 1’ (^j) 



The photon rate is the sum of (7) and (10) after insertion of k 2 = 2m 2 = 0.22g 2 T 2 (\ g + 

A»/2). 

Our result for the non-equilibrium photon rate differs from the one found in refs. [7,8], 
where the equilibrium rate was simply multiplied by the fugacities of the incoming partons 
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and the equilibrium quark mass was used as infrared cut-off. Our formula also differs from 
the one of ref. [9], where more elaborate non-equilibrium distributions (Jiittner distributions) 
were assumed, which are, however, not in line with the parton chemistry model presented 
in ref. [6]. 


IV. PHOTON SPECTRUM 


In the simplest of all models of a central ultrarelativistic heavy ion collision, one imagines 
two flat nuclei which penetrate each other and fly apart afterwards, creating a longitudinally 
expanding, cylindrical quark-gluon plasma tube between them. A Minkowski diagram of this 
event is shown in fig.2. 

For the total photon yield, we take the rate as given in the last section and the fugacities 
shown in fig.l, and integrate over the plasma space-time volume, using 


c f 


/ plasma 


T\ Dnuc 



T 0 ~Unuc 


where the times after the maximum overlap of the nuclei are r 0 = 0.7 fm/c and T\ = 4 fm/c 
and the rapidity of the nuclei is y nuc = 6 at RHIC and r 0 = 0.5 fm/c, T\ = 6.25 fm/c, and 
Unuc = 8.8 at LHC, respectively. Q is the transverse cross section of the nuclei; for gold, 
Q ~ 180/m 2 . We obtain for the photon spectra 

dn loc - rest ' 


/ 2 dn \ 
\d 2 p_\_dy j 


<y,pi 


T\ Dnuc / 

= Q j dr t J dy' \ 2E 

TO —Unuc 


d 3 pd 4 x 


( 11 ) 


I E loc ' rest _ pL cos h( 2/ - </), 
T = T(t), A = A(t) 


(The two expressions in large round brackets are Lorentz scalars; we can Lorentz transform 
them by simply transforming their argument, i.e., the photon energy. The photon energy 
in the laboratory frame equals pj_cosh(j/); the photon energy in the comoving frame of 
the plasma is then p± cosh(j/ — y r ). The superscript loc.rest serves to remind us of this 
transformation.) Here we assumed that we can use the fugacities for times r > r 0 not only 
for the central region and neglected the transverse expansion of the fireball, which can be 
treated in a hydrodynamical model [21]. 


8 



Fig.3 shows the individual contributions of different time intervals between r 0 and T\ to 
the photon spectrum. Summing up these contributions gives rise to the concave shape of the 
photon spectrum in fig.3. Clearly the spectrum is dominated by early times corresponding to 
high temperatures [8], especially for large photon energies. Hence it is rather insensitive to 
the uncertainties of the equilibration rates. Compared to chemical equilibrium X g = X q = 1 
at an initial temperature of T 0 = 300 MeV, as it was considered for SPS energies [21], 
the photon yield from the QGP is suppressed by about a factor of 10 -2 (RHIC) to 10 -1 
(LHC). This would imply that photons from the QGP cannot be observed. Although the 
photon rate is enhanced by the high initial temperature T 0 = 500 - 800 MeV, this increase 
is overcompensated by the small fugacities of our model. Using the fugacities as plotted in 
ref. [12], the photon yield is the same as in fig.3 within a factor of 2. 

Spectra for different rapidity regions are plotted in fig. 4. Within the limits given by 
the nuclear rapidities, the transverse momentum spectra do not depend on the photon 
rapidity. In other words, the photon rapidity distribution closely resembles the underlying 
quark rapidity distribution [22], This is easily understood, as the rapidity distribution of 
massless isotropic radiation is always a bell-shaped curve with a width at half maximum 
(FWHM) Ay ~ 1.6. This is small compared with the nuclear rapidities (y nuc = 6 — 8.8) 
at high energies. Therefore, if we Lorentz transform (shift) the rapidity distribution of the 
photons according to the velocity of the plasma cell they originate from, we get essentially 
the plateau-shaped distribution of the quarks, smeared out a little bit. This is schematically 
shown in fig. 5. 


V. CONCLUSIONS 

We have studied the photon production from a chemically non-equilibrated parton gas 
presumably produced at RHIC and LHC. We have neglected the pre-thermal stage as well 
as the mixed and hadronic phase of the fireball. Using the undersaturated parton densities, 
following from rate equations of the parton chemistry based on initial conditions predicted 


9 



by HIJING, the photon yield is suppressed by a factor 10 -2 (RHIC) to 10 -1 (LHC) compared 
to a fully equilibrated QGP at an initial temperature of T 0 = 300 MeV. This large difference 
relies on the fact that the HIJING model predicts a very small quark density (only 2 - 3% of 
the equilibrium value) at the time r 0 , from which on the parton distributions look thermal, 
and that the photon emission is dominated by early times. This result suggests that the 
photon emission from the plasma phase is not observable. In contrast, the Parton Cascade 
Model predicts almost complete saturation of the parton phase space densities at the onset 
of thermal equilibrium with similar initial temperatures (T 0 = 500 - 800 MeV). Hence the 
photon yield from the plasma is enhanced by about a factor of 10 2 (RHIC) to 10 3 (LHC) 
compared with our estimate and might be visible, in particular for photon energies between 
2-3 GeV [1], Consequently, photon data at RHIC and LHC might allow a distinction 
between the initial condition predicted from HIJING and the one from the Parton Cascade 
Model. 
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FIGURES 



FIG. 1. Temperature, quark fugacity, and gluon fu.gacity for an Au + Au-collision in RHIC 
(s 1 ! 2 = 200 A ■ GeV) and LHC (sf! 2 = 6000 A ■ GeV). The data was taken from [6]. 
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FIG. 2. Minkowski diagram of a relativistic heavy ion collision (RHIC) creating a quark-gluon 
plasma. We compute the photon yield of the plasma between times tq and T\. 
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FIG. 3. Midrapidity (y=0) photon yield from various stages (proper time slices) of the plasma 
e volution. The top curve is identical to the y=0 curve in the next figure. 
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FIG. 5. Rapidity distributions of photons, partons, and those photons originating from a part.on 
cell within the rapidity interval Ay part . 
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